Search for: All records

Abstract The unusual properties of superconductivity in magic-angle twisted bilayer graphene (MATBG) have sparked considerable research interest^1–13. However, despite the dedication of intensive experimental efforts and the proposal of several possible pairing mechanisms^14–24, the origin of its superconductivity remains elusive. Here, by utilizing angle-resolved photoemission spectroscopy with micrometre spatial resolution, we reveal flat-band replicas in superconducting MATBG, where MATBG is unaligned with its hexagonal boron nitride substrate¹¹. These replicas show uniform energy spacing, approximately 150 ± 15 meV apart, indicative of strong electron–boson coupling. Strikingly, these replicas are absent in non-superconducting twisted bilayer graphene (TBG) systems, either when MATBG is aligned to hexagonal boron nitride or when TBG deviates from the magic angle. Calculations suggest that the formation of these flat-band replicas in superconducting MATBG are attributed to the strong coupling between flat-band electrons and an optical phonon mode at the graphene K point, facilitated by intervalley scattering. These findings, although they do not necessarily put electron–phonon coupling as the main driving force for the superconductivity in MATBG, unravel the electronic structure inherent in superconducting MATBG, thereby providing crucial information for understanding the unusual electronic landscape from which its superconductivity is derived. 

Abstract Interactions between electrons and phonons play a crucial role in quantum materials. Yet, there is no universal method that would simultaneously accurately account for strong electron-phonon interactions and electronic correlations. By combining methods of the variational quantum eigensolver and the variational non-Gaussian solver, we develop a hybrid quantum-classical algorithm suitable for this type of correlated systems. This hybrid method tackles systems with arbitrarily strong electron-phonon coupling without increasing the number of required qubits and quantum gates, as compared to purely electronic models. We benchmark our method by applying it to the paradigmatic Hubbard-Holstein model at half filling, and show that it correctly captures the competition between charge density wave and antiferromagnetic phases, quantitatively consistent with exact diagonalization.